Precision Measurement of Cosmic-Ray Antiproton Spectrum 
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The energy spectrum of cosmic-ray antiprotons (p's) has been measured in the range 0.18 to 3.56 
GeV, based on 458 p's collected by BESS in recent solar-minimum period. We have detected for 
the first time a distinctive peak at 2 GeV of p's originating from cosmic-ray interactions with the 
interstellar gas. The peak spectrum is reproduced by theoretical calculations, implying that the 
propagation models are basically correct and that different cosmic-ray species undergo a univer- 
sal propagation. Future BESS flights toward the solar maximum will help us to study the solar 
modulation and the propagation in detail and to search for primary p components. 

PACS numbers: 98.70.Sa, 95.85.Ry 



The origin of cosmic-ray antiprotons (p's) has attracted 
much attention since their observation was first reported 
by Golden et al. H. Cosmic-ray p's should certainly 
be produced by the interaction of Galactic high-energy 
cosmic-rays with the interstellar medium. The energy 
spectrum of these "secondary" p's is expected to show a 
characteristic peak around 2 GeV, with sharp decreases 
of the flux below and above the peak, a generic feature 
which reflects the kinematics of p production. The sec- 
ondary p's offer a unique probe [0 of cosmic-ray propaga- 
tion and of solar modulation. As other possible sources of 
cosmic-ray p's, one can conceive novel processes, such as 
annihilation of neutralino dark matter or evaporation of 
primordial black holes || . The p's from these "primary" 
sources, if they exist, are expected to be prominent at 
low energies [Q and to exhibit large solar modulations 
|5|. Thus they are distinguishable in principle from the 
secondary p component. 

The detection of the secondary peak and the search 
for a possible low-energy primary p component have been 
difficult to achieve, because of huge backgrounds and the 
extremely small flux especially at low energies. The first 
and subsequent g evidence for cosmic-ray p's were 
reported at relatively high energies, where it was not pos- 
sible to positively identify the p's with a mass measure- 
ment. The first "mass-identified" and thus unambiguous 
detection of cosmic-ray p's was performed by BESS '93 
J7J in the low-energy region (4 events at 0.3 to 0.5 GeV), 
which was followed by IMAX § and CAPRICE de- 
tections. The BESS '95 measured the spectrum (To) at 
solar minimum, based on 43 p's over the range 0.18 to 
1.4 GeV. We report here a new high-statistics measure- 



ment of the p spectrum based on 458 events in the energy 
range from 0.18 to 3.56 GeV. 

Fig.p] shows a schematic view of BESS. It was designed 
pd[ and constructed Q as a high-resolution spectrom- 
eter to perform searches for rare cosmic-rays, as well as 
various precision measurements. A uniform field of 1 
Tesla is produced by a thin (4 g/cm 2 ) superconduct- 
ing coil Q, through which particles can pass without 
too many interactions. The magnetic-field region is filled 
with the tracking volume. This geometry results in an 
acceptance of 0.3 m 2 sr, which is an order of magnitude 
larger than those of previous cosmic-ray spectrometers. 
The tracking is performed by fitting up to 28 hit-points 
in the drift chambers, resulting in a magnetic-rigidity 
(R) resolution of 0.5 % at 1 GV/c. The upper and 
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FIG. 1. Cross-sectional view of the BESS spectrometer in 
its 1997 configuration. Overlayed is one of the p events. 



lower scintillator-hodoscopes provide two dE/dx mea- 
surements and the time-of-flight (TOF) of particles. The 
dE/dx in the drift chamber gas is obtained as a truncated 
mean of the integrated charges of the hit-pulses. For 
the '97 flight, the hodoscopes were placed at the outer- 
most radii, and the timing resolution of each counter was 
improved to 50 psec rms, resulting in /3 _1 resolution of 
0.008, where (3 is defined as particle velocity JT3] divided 
by the speed of the light. Furthermore, a Cherenkov 
counter with a silica-aerogel (n — 1.032) radiator was 
newly installed pa], in order to veto e~//i - backgrounds 
which gave large Cherenkov light outputs corresponding 
to 14.7 mean photo-electrons when crossing the aerogel. 

The 1997 BESS balloon flight was carried out on July 
27, from Lynn Lake, Canada. The scientific data were 
taken for 57,032 sec of live time at altitudes ranging from 
38 to 35 km (an average residual air of 5.3 g/cm 2 ) and 
cut-off rigidity ranging from 0.3 to 0.5 GV/c. The first- 
level trigger was provided by a coincidence between the 
top and the bottom scintillators, with the threshold set 
at 1/3 of the pulse height from minimum ionizing par- 
ticles. The second-level trigger, which utilized the hit- 
patterns of the hodoscopes and the inner drift chambers 
(IDC), first rejected unambiguous null- and multi-track 
events and made a rough rigidity-determination to select 
negatively-charged particles predominantly. In addition, 
one of every 60 first-level triggers was recorded, in order 
to build a sample of unbiased triggers. 

The off-line analysis pfl| selects events with a single 
track fully contained in the fiducial region of the track- 
ing volume with acceptable track qualities. The three 
dE/dx measurements are loosely required as function of 
R to be compatible with proton or p. The combined 
efficiency of these off-line selections is 83 - 88 % for R 
from 0.5 to 4 GV/c. These simple and highly-efficient 
selections are sufficient for a very clean detection of p's 
in the low-velocity (f3 < 0.9) region. At higher- velocities, 
the e~~ j n~~ background starts to contaminate the p band, 
where we require the Cherenkov veto; i.e., 1) the particle 
trajectory to cross the fiducial volume of the aerogel, and 
2) the Cherenkov output to be less than 0.09 of the mean 
output from e~. This cut reduces the acceptance by 20 
%, but rejects e _ //z~ backgrounds by a factor of 6000, 
while keeping 93 % efficiency for protons and p's which 
cross the aerogel with rigidity below the threshold (3.8 
GV/c). Fig.g] shows the /3 _1 versus R plot for the sur- 
viving events. We see a clean narrow band of 415 p's at 
the exact mirror position of the protons. The p sample is 
thus mass-identified and background- free, as the clean- 
ness of the band demonstrates and various background 
studies show. In particular, backgrounds of albedo and 
of mis-measured positive-rigidity particles are totally ex- 
cluded by the excellent /3 _1 and i? _1 resolutions. To 
check against the "re-entrant albedo" background, we 
confirmed that the trajectories of all p's can be traced 
numerically through the Earth's geomagnetic field back 



to the outside of the geomagnetic sphere. 

We obtain the p fluxes at the top of the atmosphere 
(TOA) in the following way: The geometrical acceptance 
of the spectrometer is calculated both analytically and 
by two independent Monte Carlo methods. The live 
data-taking time is directly measured by two independent 
scaler systems gated by the "ready" gate which controls 
the first-level trigger. The efficiencies of the second- level 
trigger and of the off-line selections are determined by us- 
ing the unbiased trigger sample. The TOA energy of each 
event is calculated by tracing back the particle through 
the detector material and the air. The interaction loss 
of the p's is evaluated by applying the same selections to 
the Monte Carlo events generated by GEANt/gheisha, 
which incorporates |L6|] detailed material distribution and 
correct p-nuclei cross sections. We subtract the expected 
number |17]] of atmospheric p's, produced by the colli- 
sions of cosmic-rays in the air. The subtraction amounts 
to 9 ± 2 %, 15 ± 3 % and 19 ± 5 %, respectively, at 0.25, 
0.7 and 2 GeV, where the errors correspond to the max- 
imum difference among three recent calculations fl7j-|l9]l 
which agree to each other. Proton fluxes are obtained 
in a similar way. Atmospheric protons are subtracted by 
following Papini |2(j] . 

Table | contains the resultant BESS '97 p fluxes and 
p/p flux ratios at TOA. The first and the second errors 
represent the statistical |2^] and systematic errors, rc- 
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FIG. 2. The identification of p events. The solid lines 
define the p mass band used for the spectrum measurement. 



TABLE I. Antiproton fluxes (in xlCT 2 nrW^GeV" 1 ) and p/p ratios (in xl0~ 5 ) at TOA. T (in GeV) define the 
kinetic energy bins. Np and Tp, respectively, are the number of observed antiprotons and their mean kinetic energy in each 
bin. The eighth bin of BESS '95 flux actually covers the energy region from 1.28 to 1.40 GeV. 
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spectively. We checked that the central values of the 
fluxes are stable against various trial changes of the 
selection criteria, including uniform application of the 
Cherenkov veto also to the low (3 region. The dominant 
systematic errors at high and low energies, respectively, 
are uncertainties in the atmospheric p calculations and 
in the p interaction losses to which we attribute ±15 % 
relative error. As shown in Table Q, the BESS '97 fluxes 
are consistent [g2[ with the '95 fluxes in the overlapping 
low-energy range (0.2 to 1.4 GeV). The solar activities 
at the time of the two flights were both close to the min- 
imum as shown by world neutron monitors and by the 
low-energy proton spectra |23] measured by BESS. 

Shown in FigJ| is the combined BESS ('95+'97) spec- 
trum, in which we detect for the first time a distinctive 
peak at 2 GeV of secondary p, which clearly is the dom- 
inant component of the cosmic-ray p's. 

The measured secondary p spectrum provides crucial 
tests of models of propagation and solar modulation, 
since one has a priori knowledge of the input source 
spectrum for the secondary p, which can be calculated 
by combining the measured proton and helium spectra 
with the accelerator data J24| on the p production. The 
distinct peak structure of the p spectrum also has clear 
advantages in these tests over the monotonic (and un- 
known) source spectra of other cosmic-rays. 

The curves shown in Fig.[3| are recent theoretical cal- 
culations for the secondary p in diffusion model | J25|]26| 
and leaky box model |27[|2(j| , in which the propagation 
parameters (diffusion-coefficient or escape length) are de- 
duced by fitting various data on cosmic-ray nuclei such as 
Boron/Carbon ratio, under the assumption that the dif- 
ferent cosmic-ray species (nuclei, proton and p) undergo 
a universal propagation process. These calculations 



use as essential inputs recently measured proton spectra 
p3| , pl 32|, which are significantly (by a factor 1.4 to 1.6) 
lower than previous data in the energy range (10 to 50 
GeV) relevant to the p production. 

These calculations reproduce our spectrum at the peak 
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FIG. 3. BESS '95+'97 (solar minimum) antiproton fluxes 
at the top of the atmosphere together with previous data. 
The error-bars represent the quadratic sums of the statistical 
and systematic errors. The curves are recent calculations of 
the secondary p spectra for the solar minimum period. 
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FIG. 4. BESS '95+ '97 p/p flux ratios with previous data 
[9] and the theoretical calculations. 

region remarkably well within their ±15 % estimated ac- 
curacy p8| . This implies that the propagation models 
p3[ are basically correct and that different cosmic-ray 
species undergo a universal propagation process. 

At low energies, the calculations predict somewhat di- 
verse spectra reflecting various uncertainties [ B4J , which 
presently make it difficult to draw any conclusion on pos- 
sible admixture of primary p component. As noted in 
Ref. ||, the rapid increase of solar activity toward the 
year 2001 will drastically suppress the primary ^compo- 
nent such as from the primordial black holes |3J, while 
changing the shape of the secondary p only modestly 
p7| . This will help us to separate out the "primary" 
and "secondary" components at low energies. Future an- 
nual BESS flights will thus be important to search for 
primary p component, to study solar modulation, and to 
investigate further details of the propagation. 

Since most previous data were presented in the form 
of p/p flux ratios B, a compilation is made in Figfl, 
which shows again the unprecedented accuracy of our 
measurement. 
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